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ON NUMERICAL DENSITY APPROXIMATIONS OF 
SOLUTIONS OF SDES WITH UNBOUNDED COEFFICIENTS 

LINGHUA CHEN, ESPEN ROBSTAD JAKOBSEN, AND ARVID NAESS 


Abstract. We study a numerical method to compute probability density 
functions of solutions of stochastic differential equations. The method is some¬ 
times called the numerical path integration method and has been shown to be 
fast and accurate in application oriented fields. In this paper we provide a rig¬ 
orous analysis of the method that covers systems of equations with unbounded 
coefficients. Working in a natural space for densities, , we obtain stability, 
consistency, and new convergence results for the method, new well-posedness 
and semigroup generation results for the related Fokker-Planck-Kolmogorov 
equation, and a new and rigorous connection to the corresponding probabil¬ 
ity density functions for both the approximate and the exact problems. To 
prove the results we combine semigroup and PDE arguments in a new way 
that should be of independent interest. 


1. Introduction 


Over the past decades there has been a large number of publications in the field of 
stochastic dynamics and its various application areas - including physics, biology, 
engineering, and finance [34, 2, 11]. In this field, the response of dynamical systems 
to stochastic excitation is studied, and the typical model is a (system of) stochastic 
differential equations (SDEs) of the form 



where b : K+ x —>• a : R+ x R'^ —R'^^”, Bt is an n-dimensional Brownian 

motion, and the initial data Z is a random variable in R*^ independent of Bt- The 
solution Yt of the SDE is a state space process in R^. 

In most cases the solutions of such problems must be computed numerically, and 
various discrete approximation methods are widely used in many application areas 
[24], There are two main approaches: (i) Path-wise approximations of the SDE 
based on stochastic simulation, and (ii) approximations of the statistics or distri¬ 
butions of the SDE. The first approach is more efficient in high dimensions and the 
second in low dimensions. For path-wise approximations we refer to [24] for the 
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classical literature and e.g. to [31] and references therein for some promising re¬ 
cent developments. Approach (ii) is deterministic and based on approximating the 
forward or backward Kolmogorov partial differential equations. In this paper we 
study a method of the second type, called the numerical path integration method 
[36, 29, 39, 41, 40, 10] or density tracking method [5]. This method is an explicit 
deterministic iteration scheme that produces approximate probability density func¬ 
tions (PDFs) for the solution of the SDE. The iteration step is based on a short time 
approximation of the SDE. For simplicity we use here the strong Euler-Maruyama 
method, the most basic numerical method for SDEs. The convergence of our path 
integration method is therefore equivalent to the convergence of the PDEs of the 
Euler-Maruyama method. 

The path integration/density tracking approach (i.e. simulating the PDFs) enjoys 
several favorable properties. First, it introduces an extra perspective to the system, 
which enables deeper insights and invites broader mathematical tools. Secondly, as 
an explicit method, one can formally implement the path integration algorithm on a 
vast number of scenarios. Since the formulation is deterministic, it is also free from 
perturbation by extreme outcomes during stochastic simulation. Finally, the result 
of the method is an explicit density function rather than bundle of random paths. 
This means that many characteristics of the system become more transparent, and 
can be captured and displayed by e.g. visualisation methods. 

The numerical path integration method has been applied in many fields, including 
financial mathematics [29, 39, 40, 10]. Many of these studies show that it can 
provide highly accurate numerical solutions [36, 41, 35, 5]. Even though convergence 
problems have been reported in some cases, cf. Section 7.1 in [35], little or no 
emphasis has been devoted to conditions for convergence of the method in the path 
integration literature. A very natural and relevant mode of convergence for this 
method is strong L^-convergence of the resulting densities. Such convergence seems 
not to be a direct consequence of either strong or weak convergence of the Euler- 
Maruyama scheme. In fact, -convergence of densities implies weak convergence of 
processes but the converse is not true in general. Nevertheless, very recently there 
are some results for one-dimensional problems in [5], and when the coefficients are 
bounded, L^-convergence results follow from some of very precise error expansions 
for densities of the Euler-Maruyama method, e.g. in [4, 20, 19], or from so-called 
Feynman formulas in mathematical physics [8]. More on this below. However 
from an application point of view, it is important to consider also models with 
unbounded coefficients, since they naturally appear in many of the papers and 
applications mentioned above. ^ 

The aim of this paper is to provide a rigorous analysis of method in the path inte¬ 
gration setting along with the T^-convergence of the PDFs it produces in the case 
of unbounded coefficients. Specifically, we obtain stability, consistency, and new 
convergence results for the method, new well-posedness and semigroup generation 
results for the associated Fokker-Planck-Kolmogorov (FPK) equation, and a new 
and rigorous connection to the corresponding probability density functions for both 


^By a change of variables, it is often possible to reduce to a problem with bounded coefficients 
and then use existing convergence results. However going back to the original variables, we 
would then get a weaker form of convergence, one that no longer implies convergence. 
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the approximate and exact problems. In particular, our semigroup generation 
result for our rather general elliptic FKP operator seems to be new. 

We also investigate in details the method applied to the Ornstein-Uhlenbeck process 
which has unbounded drift. Using explicit transition densities and direct compu¬ 
tations, we show that the method (i) converge in the general case and (ii) converge 
with a linear rate if additional moment conditions are assumed. The second result 
is consistent with results for bounded coefficients in e.g. [4, 20, 19]. Even though 
we do not have a counter example, the very explicit proof seems to indicate that 
convergence rates can only be obtained under some additional moment assumptions 
on the densities. 

We study the numerical path integration method from a PDE and semigroup point 
of view, an approach which seems to be new to this setting at least when it comes 
to the convergence results. The main idea is to show that (i) the EPK operator 
corresponding to the SDE generates a semigroup, (ii) the iteration of the discrete 
path integration operator does converge to this semigroup (the solution of the EPK 
equation), and finally (iii) that convergence of the densities follows from these 
results. For (i), we prove that the FKP operator and its adjoint are dissipative and 
then use the Lumer-Phillips theorem to prove that the FPK operator generates a 
strongly continuous contraction semigroup on L^(IR‘^). To prove the convergence 
results in (ii), we use the Chernoff product formula which is a generalisation of the 
well-known Kato-Trotter product formula. Note that our approach to (ii) is similar 
to the approach for Feynman formulas in [8]. 

A main difficulty is then to prove the semigroup generation in our setting of 
unbounded coefficients. Classical results usually assume the coefficients to be 
bounded, cf. the discussion in Chapter 11 of [28]. Further complications comes 
from the fact that we have to work in the space which is not reflexive and the 
adjoint space L°° which is not separable. In fact one of the most difficult parts of 
our paper is to show that the adjoint operator is dissipative on L°°. To do this we 
develop a new and non-trivial argument using ideas from viscosity solution theory 

[12] , the weak Bony maximum principle [6, 30], and recent elliptic regularity results 
from [42]. This argument could be of independent interest. 

Other authors have obtained semigroup generation results for similar problems 
with different techniques. Yet, most work seems to be devoted to the generation 
problems for the adjoint of the FPK operator in various situations. For instance 

[13] is concerned with the space of bounded (continuous) functions and unbounded 
coefficients, while in [1] the authors give results in for bounded coefficients. The 
author in [9] considered degenerate operators where the coefficients have bounded 
first and second derivatives. The result closest to our case seems to be found in [15] 
where the drift term had to be dominated by the ’’square root” of the diffusion term 
and no zero order term is allowed. A semigroup generation result of the adjoint 
operators on is proved there. However, we cannot use the results in [15] in this 
paper since our non-divergence form FPK operators are equivalent to divergence 
form operators with additional unbounded zero order terms. 

In addition to providing rigorous results for the numerical path integration method, 
this paper seems to represent the first attempt to use semigroup and PDE arguments 
to obtain strong convergence results for probability densities of SDE approximations 
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(but see also [8]). It is an interesting question how far such methods can be pushed 
compared to other methods. As far as we know, two other methods are described 
in the literature: A parametrix approach of Konakov et al. [25, 26, 33], and an 
approach based on Malliavin calculus developed by Talay, Bally, et ah, see e.g. [4, 
20, 3] and references therein. In the parametrix approach, a weighted L°“-estimate 
on the convergence of (transition probability) densities, as well as non-asymptotic 
bounds for the densities [33], is obtained under uniform ellipticity assumption and 
boundedness of coefhcients. In the Malliavin approach various error expansions and 
estimates are obtained under boundedness, certain smoothness, and ellipticity or 
Hormander conditions. 

Under additional assumptions, the -convergence of densities of the Euler-Maruyama 
scheme can be obtained from some of these Malliavin results. In [4] a linear conver¬ 
gence rate was obtained for SDEs with smooth and bounded coefficients^, but the 
behavior of the estimate is unclear for small time. In [20] the estimate was improved 
but it would still blow up as t 0. Later the authors of [19] derived very accurate 
estimates that would lead to a linear rate which is uniform for small time. At 
the same time they lowered the regularity requirements of the coefficients to C^’^. 
For the Euler-Maruyama method, the estimates of [4, 20, 19] requires bounded co¬ 
efficients and uniform ellipticity. According to Remark 2.22 of [3], the ellipticity 
requirement can be relaxed to a so-called weak Hormander condition when coeffi¬ 
cients belong to . In this case, the results of [3] give -convergence for fixed 
times (but no estimates) of densities/PDFs for the Euler-Maruyama scheme. 

To summarise, in our approach to study properties and convergence of the numerical 
path integration method we obtain the following results: New semigroup generation 
and well-posedness results for the FPK equation, allowing for coefficients with linear 
growth and differentiability, no invariant measure is required, the growth of drift 
and diffusion terms are not necessarily related to each another, and the initial data 
is not required to have any finite moments or differentiability. Under the same 
conditions, we prove convergence, locally uniformly in time (so no deterioration 
for small time). To our knowledge, this seems to be the first general strong L^- 
convergence result for densities/PDFs of SDEs with unbounded coefficients. Even 
though we do not obtain any rate of convergence or error expansion in general, by 
studying the Ornstein-Uhlenbeck process, we find indications that no error estimate 
may exist without additional (moment) assumptions on the densities/PDFs. Hence 
when the coefficients are unbounded, the best result for general densities could be 
mere convergence without any rate. 

We also emphasise that our results connects the mild solution of the FPK equation 
and the PDF of the SDE, and that existence and convergence of the latter follows 
from the existence and convergence (by semigroup methods) of the former. More 
refined existence results for PDFs have been obtained by probabilistic methods, we 
refer to e.g. [7, 16, 21, 22] and references therein. 

Layout and Notation. The rest of this paper is laid out as follows. In Section 
2 we state our assumptions and the main results. Included is a discussion of the 
connection between SDEs, FPK equations and densities and the definition of the 

^The coefficients are bounded in [4]. This follows from their condition (H), estimate ( 6 ), and 
footnote 1. In particular, in [4] the number 70 , introduced by (1.10) in [27], is zero. 
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path integration method. At the end of the section, we discuss our results and give 
some examples. The rest of the paper is then devoted to the proofs of these results. 
In Section 3 we prove the well-posedness of the FPK equation, the connection 
to probability densities, and the convergence result - assuming the semi-group 
generation and a strong T^-consistency result for the path integration method, 
which are proved in Section 4 and Section 5, respectively. At the end. Section B is 
an appendix to the calculation of the convergence rate of Example 2.13. 

The following notation and abbreviations are used H-H^^ := IHloo := 

||•|lioo(Rd), dt := V := =: (^iv ,ddY, essinf is the essen¬ 

tial infimum, E denotes the expectation; C^, C^, V the spaces of functions with 
bounded continuous derivatives up to A:—th order, smooth compactly supported 
functions, and distributions {V is the dual of C^) respectively, (^([O, T]; L^(]R‘^)) 
the space of functions u{x,t) such that sup^gjoj^j ||M(-,t)||i < oo, and for all t € 
[0, T], ||it(-, s) — u(-, t)||i —0 as s —>■ t. PDF - probability density function, SDE - 
stochastic differential equation, FPK - Fokker-Planck-Kolmogorov. 


2. Main Results 

In this section we formulate the discrete path integration method, state the as¬ 
sumptions and main results, and provide a discussion and examples. We restrict 
ourselves to time-homogeneous SDEs, i.e. equations of the form 


( 2 . 1 ) 


dYt = b{Yt) dt + a(Yt) dBt, 

Yo = Z, 


where b = [bi, - ■ ■ , bj)"^ : —>■ and a = {<Jij)dxn '■ t are functions, 

Bt a d-dimensional Brownian motion, and Z is a random variable in R'^. 

Under suitable assumptions (cf. [37] or [38]), the solution Yt of (2.1), is a Markov 
process, with infinitesimal generator A* defined as 

d 1 ^ 

(2.2) := ^^bi{x)dicj){x) + - ^ aij{x)didj(j){x) for all (j) G C'))°(R'^), 


Li=i 


where 


a :— (jfj — {ciij'jdxd- 

Moreover, Yt has a PDF u{t,x) - a non-negative function u such that 

(2.3) ^(t){Yt) = [ <j){x)u{t,x)dx for all (j) G 

Jk‘‘ 

Via an adjoint argument, u formally satisfies the FPK equation 
dtu{t,x) = Au{t,x), X G R'^, t > 0, 

u{0,x) = uo{x), a; G R'^, 

where Uq is the PDF of Z and A the adjoint of A*, 

d 1 ^ 

Av{x) := -'^di{bi{x)v{x)) + 2 X! didj{aij{x)v{x)). 

i—1 i,j — l 


(2.4) 


(2.5) 
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The evolution of the u can be regarded as the action of the propagation operator 
(or semigroup) Ps defined by 

(2.6) u{y,t + s) = [Psu{-,t)]{y) for all t,s>0, ?/G 

Let us state our assumptions. 

(Cl) 6 : R'^ —R'^ and a : R'^ —R'^^" are functions, and there exists a 
constant K > 0 such that for all x G R^, j = 1,... ,n, and i,k = 1,..., d, 

\dk<Jtj{x)\ + \dkbi{x)\ < K. 

(C2) There exists a > 0 such that for all x, y G R^, 

y^a{x)y > a|yp. 

(C3) Z and Bt are independent, and Z has a PDF uq, i.e. uq > 0 and ||mo|Ii = 1- 

Remark 2.1. Standard strong well-posedness of the SDE (see below) requires only 
Lipschitz coefficients. The additional regularity and the uniform ellipticity condi¬ 
tion (C2) is needed to study the PDFs of the SDEs and path integration method 
(cf. proofs of Proposition 2.10 and Lemma 5.8). Similar but stronger assumptions, 
including bounded coefficients, are used by all other papers discussing densities of 
SDE approximations, see the introduction and e.g. [19]. 

The following result is then classical, cf. Theorem V.7 in [38]. 

Proposition 2.2. Assume (Cl) and (C3). Then the SDE (2.1) has a unique strong 
solution Yt- 

Now we give our first main result - well-posedness for the FPK equation (2.4). 

Theorem 2.3. Assume (Cl), (C2), and uq G L^(R'^). Then the FPK equation 
(2.4) has a unique (mild) solution u G (^([O,Tj;L^(R'^)). 

The derivation of FPK equation (2.4) was formal. Yet our second main result 
confirms that the probability distribution of (2.1) coincides with the (mild) solution 
of FPK equation (2.4). 

Theorem 2.4. Assume (Cl) - (C3), and letYt andu{t,x) be the solutions of (2.1) 
and (2.4) respeetively. Then (2.3) holds and u{t,-) is the PDF ofYt. 

The proofs of Theorems 2.3 and 2.4 will be given in the next section. 

Remark 2.5. (a) The proof shows that Theorem 2.3 still holds if the regularity 
in assumption (Cl) is reduced to 6 G and a G C^. 

(b) A difference between Theorem 2.3 and existing results is that we allow (and 
need to allow) the zero order coefficient to be unbounded, cf. e.g. Section 4 in [15]. 


(c) The solution u we obtain in Theorem 2.3 is a mild solution of (2.4) (cf. [14]). 
In fact we have more regularity than stated since u(-,t) belongs to the domain of 
A. A limit argument immediately shows that u is also a distributional solution of 
(2.4). In view of uniform ellipticity of A and smoothness of its coefficients, it is 
standard to prove higher regularity and that u is a classical solution. 
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We now introduce and analyse the discrete path integration method. It is a method 
for computing approximations of the PDF of solution process of (2.1) by: 

(i) Finding a short-time approximation Pr of the propagator/semigroup Pt, and 


{a) computing the approximation at time t = nr using the explicit iteration 

(2.7) Un+l{x) = {PrUn){x) for 77 = 0,1,.... 

The most obvious way to find a Pt is through a time-discretization of the underlying 
SDE. For simplicity, we consider here the Euler-Maruyama method [24], the most 
basic and widely used SDE approximation scheme: Fix a time step r > 0 and let 
ABn := Bnr — and define the Markov chain approximation Xn of the 

solution process Yt of the SDE (2.1) by 

{Xo = Z, 

( 2 . 8 ) 

I Xn+l — Xn + b{Xn)T -f a{Xn)ABn, 77 — 0, 1, ... . 


In the rest of the paper we will use the following notation. 

Definition 2.6. (i) Un is the PDF of the solution Xn of (2.8). 

(m) Pr is the one step propagator of i.e. 

Un+i{y) = {PTUn){y) for all y e r > 0, n > 0. 


(in) k{y, x, r) is the transition kernel of P,- ” a non-negative function satisfying 
(2.9) iPTUn){y) = / k{y,x,T)un{x)dx for all y G r > 0, n > 0. 

Since AP„ is a multivariate Gaussian variable, we have the following classical result. 

Lemma 2.7. Assume a := cra'^ is strictly positive definite, and let k be given by 
Definition 2.6. Then for all x,y G and r > 0, 

Pyk^^, X, t) = exp I i^^{x + b{x)T) - a{x)^ 


where T denote the Fourier transform, and 

exp{-5(y- {x + b{x)T))'^ (ra)-^ {x){y - {x + b{x)T))} 


(2.10) k{y,x,T) = 


(27r)^/'^| det(Ta(a;))|^/^ 


Existence and uniqueness of the solution Xn come for free, since the iteration (2.8) 
is explicit. In view of Lemma 2.7, we then have the following result. 

Proposition 2.8 (Well-posedness). Assume (Cl) - (C3). Then there exists a 
unique PDF Un of the solution Xn of (2.8). Moreover, Un is explicitly given by 
(2.7), (2.9), and (2.10). 


Note that (C2) implies that a{x) is strictly positive definite. A classical computation 
using (2.10) then shows that ||fc(-,a:,r)||i = 1, and hence (see the next Proposition) 
that ||Pr|| = 1. Proofs can be found e.g. in the discussion of Markov operators in 
[28]. 
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Proposition 2.9 (L^—stability). Assume (C2), t > 0, u G and Pt is 

given by Definition 2.6. Then ||Pt-m|1i < ||u||i. 


By linearity, we immediately get continuous dependence on the initial data, 

and by iteration that the discrete path integration method is L^-stable: 

||un||i < ||wo||i for all n = l,2,.... 


Next we study the L^-consistency of the discrete scheme. 


Proposition 2.10 (L^— consistency). Assume (Cl), (C2), r > 0, andu G 
Then there exists a constant C > 0 such that 



— Au 


1 


< Ct. 


We prove this result in Section 4. Higher regularity of 6 or ct will not improve 
the (linear) rate in Proposition 2.10. It is the maximal rate for the problem, cf. 
Proposition 2.14 and its proof. 

Finally we state the third main result of this paper, the -convergence of the 
discrete path integration method. 

Theorem 2.11 (L^-convergence). Assume (Cl), (C2), uq G L^(IR‘^), u{t,x) is the 
solution of FPK equation (2.4), and define Un(t,x) := {Pf)^uo){x). Then 

(2.11) lim sup |lM(-,t) - M„(-,t)||i = 0. 

n^oo 

Remark 2.12. (a) Since the PDFs in general are (non-negative) L^-functions, strong 
L^-convergence is a very natural mode of convergence to consider. 

(b) By assumption the initial PDF ug is only required to belong to L^. It is not 
required to have any finite moments or differentiability to obtain an -convergence, 
and the convergence is uniform in time on [0,T] for any T > 0. 

(c) -convergence of PDFs is not a direct consequence of either strong or weak 
convergence of the solution process. It is strictly stronger than weak convergence 
of the corresponding process by Proposition A.l and Example A.2 in Appendix A. 


When the coefficients are also bounded, many authors have obtained not only 
convergence but even error estimates and error expansions for densities of the Euler- 
Maruyama scheme (see the introduction). 

In the unbounded coefficients case however, there may not be any general error 
bound without extra moment assumptions on the density, and hence mere conver¬ 
gence would be the best one can hope for. A first indication of this appears already 
in the Gaussian bounds on the derivatives of the SDE’s transition probability den¬ 
sity p{t,x,y) in [27]. When the coefficients are unbounded, these bounds are no 
longer bounded in x, and hence can no longer be used to derive the same estimates 
of approximations of p{t,x,y) as in the bounded case [4, 19]. Another indication 
is given by Example 2.13 below. In this example we study the Ornstein-Uhlenbeck 
process. We find both the exact and approximate transition kernels and compute 
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explicitly the conditions for having either mere convergence in or a linear rate 
for the corresponding PDFs. In the latter case, we find that additional moment 
assumptions are needed (cf. Proposition 2.14). 

Example 2.13. The Ornstein-Uhlenbeck process Yt is given by 

(2.12) dYt = bYtdt + adBt and Yq = x, 

where b, a gM. and 1), Bt are processes in R. The unique strong solution is 


Tt = 


a e 
Jo 


— bs 


dB, 


Due to the non-stochastic integrand in the Ito-integral it is easy to see that this 
process has a Gaussian law (cf. e.g. Section 11.5 in [14]), 

^2 


N (xe“, ^ (e» - 1)) , 


and hence the transition kernel is given by 

1 


(2.13) 


kiy,x,t) = 


2^|l (e2bt _ 1) 


: exp < - 


2^ (e2^* - 1) 


The corresponding Euler-Maruyama scheme with t = t/n is 

\Xo = x, 

= Xn + bXnT + CrAi?„, 
and has the following one step transition kernel by Lemma 2.7, 


(2.14) 


(2.15) 


Ky,x,T) = 


1 




: exp 


(y - (1 + bT)xf 
2rcr2 


For t = T small, it is clear that k[y, x, r) is close to k(y, x, t). 

For the convergence, we have the following result. 

Proposition 2.14. Assume 5, tr G R, T > 0, and uq G L^(R). Then 


lim sup 

ri-)-oo tg[Q T'] 

If in addition, either 

/ \xuo{x)\dx < oo 


Ptuo - Pf/^Uo 


= 0 . 


\xuQ{x)\dx < oo. 


then there exist C, N > 0, only depending on {b,a,T,uo), such that 


PtUo P^^^Uq 


Ct 

< — for all n > N. 
1 n 


The first part of Proposition 2.14 is consistent with Theorem 2.11, while the sec¬ 
ond part shows that a linear rate can be obtained under additional assumptions. 
The rate is consistent with estimates for SDEs with bounded coefficients, see e.g. 
Corollary 2.7 in [4] and Theorem 2.3 in [19]. 

We give a direct proof of Proposition 2.14 in Appendix B. Even though we do not 
have a counter example, the very explicit proof seems to indicate that convergence 
rates can only be obtained under additional assumptions on the densities. 
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Remark 2.15. The Cauchy distribution uq{x) = has no moments, but still 

satisfies /j^ \xuQ{x)\dx < oo 

3. Proofs of Theorems 2.3, 2.4, and 2.11 

We first prove Theorems 2.3 and 2.11 since the latter is needed in the proof of The¬ 
orem 2.4. We will use semi-group theory and we refer to [14] for more information 
and precise definitions of the concepts used below. A crucial step in the proof is 
to obtain the new generation result, Theorem 3.1. From the generation result we 
obtain a solution of the FPK equation (2.4), and Theorem 2.3 follows. Then we 
prove that FPK operator A is dissipative. Armed with generation and dissipativity, 
we use the Chernoff product formula to show that the solution of the discrete path 
integration method converges to the solution of the FPK equation (2.4) and hence 
Theorem 2.11 follows. Finally we show Theorem 2.4, that the solution of (2.4) is 
just the PDF of the solution of (2.1). 

In the rest of this section we let A be defined in (2.5) and take D{A) := 

First we state the generation result. 

Theorem 3.1. Assume (Cl) and (C2). Then {A,D{A)) is closable and its closure 
{A,D{A)) generates a contraction semigroup Pt on (L^(]R‘^), H-jlj^). 

This result is relatively technical, and the proof is left to Section 5. 

Lemma 3.2. Assume (Cl) and (C2). Then the semigroup Pt from Theorem 3.1 
is a strongly continuous semigroup on (L^(]R‘^), H-jl]^). 

Proof. The result follows by Proposition 1.5.3 in [14] with Tt = Pt if we can verify 
that there exist S > 0, M > 1 and a dense subset D C X such that 

(1) llTtll <MforalHe [0,5], 

(2) lim Ttx = X for all x € D. 

Since Pt is a contraction semigroup by Theorem 3.1, condition (1) is satisfied with 
M = 1. Moreover, assumption (2) holds on the dense set D = D{A), since by the 
definition of the generator A: 

lim (PtX — x) = lim t (Ptx — x) /t = f lim f | lim - (PtX — x) = 0 ■ Ax = 0. 

t-i-O-l- ^ ’ i-*-o+ ^ ' \t^o+ J i-s-0+ t 

Hence (1) and (2) holds, and we conclude that Pt is strongly continuous. □ 

Proof of Theorem 2.3. In view of Theorem 3.1 and Lemma 3.2, the closure of 
{A,D{A)) generates a strongly continuous semigroup Pt on (L^(R'^), jj-jj^). Let 
u{t,x) := {PtUo){x), and note that u £ (^([O, T]; L^(]R^)) by strong continuity. 
Then by e.g. Proposition II.6.4 in [14], u is the unique mild solution of the FPK 
equation (2.4). □ 

To proceed, we introduce the following definition. 

^An absolutely continuous L^-functions not satisfying either of the conditions: uq(x) = 
Xlk (*)(1 - - fc|), where R := [k - 1/k^ ,k -|- 1/fe^] and k {0}. 
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Definition 3.3. A linear operator {B,D{B)) on a Banach space (X, H-H) is dissi¬ 
pative if 

|l(A-BH|>A|lii|| 

for all A > 0 and all u € D{B). 

One can refer to Section 11.3 of [14] for more discussion of this concept. 
Proposition 3.4. Assume (Cl). Then the operator {A,D[A)) is dissipative on 

(iHR"),IMIi)- 


Proof. Let u G D{A) and define := {a; G : ±u{x) > O}. Since C 
juj = ±u on E^, and |(A — A)u\ > (A — A)(±m), 


||(A-A)u||^ > / \{X-A)u\>f (A-A)|u|. 

Je+ue~ Je+ue~ 

A\u\ < 0, 




We claim that 
(3.1) 

and hence the proposition follows since 

||(A-A)m||i > f Ajuj = A||u||i. 

JE+UE- 

Now we prove claim (3.1) for the E^ case. If A+ ^ 0, we can approximate it by 
sets Eff^ = {x : u{x) > £„}, 0 < e„ —0, with boundaries. This can be done 
since by Sard’s theorem and the implicit function theorem, A+ has boundary 
for a.e. 0 < e: < maxM+. Note that Ef C E^ C suppit which is compact. 

Then we write A in divergence form, cf. (2.5), 

A<j) = - div [aV(j) + (div oi, • • • , div OdY'(f — 2b(j)) , 

and use the divergence theorem. We consider first the diffusion term, 

/ div (oVm) = / aVu • n, 

Je+ JdE+ 

where n the exterior unit normal vector of dE)'. Since dEf' is an e-level set of u 
and u is decreasing in the outward direction at 9A+, Vu = — fin for some /3 > 0. 
Then since a = aa'^ is positive semi-definite 

aVu ■ n = —/?n^an < 0, 

and hence 

/ div (aSJu) < 0. 

J Et 

Next we estimate the convection part. Since m = e: on dEf, 

/ div( (divoi, • • • ,divad)^M — 26u) = / ((divoi, • • • , diva^)^ — 26)u • n 

Je+ ^ ' JdEt 


= e 


/ ((div m, • • • , div Od)^ — 26) • n = e: / div( (div m, • • • , div a^)^ — 26), 
JaE+ JEt ^ ' 
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which is bounded by Ce for some C = C{a, b, suppit) since a and b are smooth and 
u has compact support. Hence 

(3.2) [ Au< Ce. 

JEt 

Note that x^+ —>■ Xe+ a-®- Since Au belongs to we may then use the 

dominated convergence theorem to pass to the limit in (3.2) as e —>■ 0^ and obtain 
(3.1). The E~ case is similar and will be omitted. The proof is complete. □ 


Now we proceed to prove Theorem 2.11. To do that we need the Chernoff Product 
Formula (e.g. Theorem 111.5.2 in [14]): 

Theorem 3.5. Let 'K be a Banach space, V : R+ —>■ /1(X) a function such that 
H(0) = id and ||H"'(t)|| < M for all t > 0, all n £ N, and some M > 1. Assume 


(3.3) 


Bu := lim 

t —^0-|- 


V(t)u — u 
t 


exists for all u £ D G where D and (A — B)D are dense subspaces in X for some 
A > 0. Then the closure B of B generates a bounded strongly continuous semigroup 
{Tft) : t > 0}, given by 

(3.4) T{f)u= lim V'^{t/n)u, 

n—¥-\-oo 


for all u £'K uniformly for t £ [0, T]. 


This result does not require the approximation operators V ft) to be a semigroup. 
This is important since the discrete path integration operators {A-} are bounded 
on the Banach space L^{W^) by Proposition 2.9. But in contrast to their continuous 
counterpart, they no longer form a semigroup of operators^. 


Proof of Theorem 2.11. We use Theorem 3.5 with X = L^(R'^), V{h) = Ph, B = A, 
and D = D{A). Let us verify the conditions. By Proposition 2.10, 

An = lim BleJi 

/i->o h 

for all u £ D{A). By the definition of D{A) and Theorems 3.4 and 3.1, A is densely 
defined, dissipative, and the closure A generates a contraction semigroup Pt on 
L^(R‘^). Then by the Lumer-Phillips Theorem (e.g. Theorem 11.3.15 in [14]), this 
implies that (A — A)D{A) is dense in L^(]R‘^) for some A > 0. 

Since the semigroup is strongly continuous by Lemma 3.2, we have verified all the 
conditions for Theorem 3.5. Hence hence can conclude that 


as n 


\WitG)-Un\\i= 

oo uniformly in [0,T]. The proof is complete. 


0 , 


□ 


'*E.g. in Example 2.13, Pf ^ P 2 t since 


{P?my) = 


{P 2 -rf)(y) = 


f{P) 


x/47TTa^ {1 + l‘£) 


exp< - 


(y — (1 + 2bT + b‘^T^)xy 

+ b^r'^12) 


J* dx, 


/W 


■s/47T 


: exp < - 


(y - (1 + 26r)x)2 


4ra^ 


dx. 
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Now we prove Theorem 2.4, using among other things Theorem 2.11. 


Proof of Theorem 2.4- We will show that (2.3) holds. 

Assume first that (j> € In view of the Euler-Maruyama scheme (2.8) cor¬ 
responding to Yt with T = ^, the PDF of its solution is := and 

hence 

(3.5) E^(A„) = j (j){x)un{x)dx for any ^ € C°(]R'^). 

dRd 

The right hand side is close to what we want, since by Theorem 2.11, 

< ll<^lloo l|w — Wnlli —>■ 0 as n —>■ oo. 

te[o.T] 



To see that the left hand side of (3.5) is also close to the sought after expression, 
we will use the weak convergence of the Euler-Maruyama method. Since our initial 
distribution uq is not assumed to have second moments, we need to introduce the 
following auxiliary SDE and Euler-Maruyama scheme for a fixed x G K: 

f dVt = b(Yt) dt + cr(Yt) dBt, f Xn+i = A„ -|- b{Xn)T a{Xn)XBn, 

|yb = X, |Ao = X. 

Both problems have unique (strong) solutions (Proposition 2.2), so we may define 
v{t,x) ■.= ¥.4>{Yt) and u„(x) := E(/)(A„). 

By the strong Markov property of solutions of SDEs, cf. Theorem V.32 in [38], 
E(^(rt) = E [E (</.(Tt)|ro = ^)] = E Kt, Z)]. 

Similarly, since {Xn} is a Markov chain, 

E0(A„) = E [E (0(A„)|Ao = Z)] = E K(Z)]. 

Therefore 

|E^(Yt) - E(^(X„)| < E |i;(t, Z) - v„(Z)| = f \v{t,x) - Vn{x)\uo{x)dx. 

jR-i 

The latter integrand is dominated by the integrable function 2 ||^||^, and by the 
weak convergence of the Euler-Maruyama scheme, cf. Theorem 14.1.5 in [24], 

\v{t, x) — Vn{x)\ —>■ 0 as n —oo for every x G 

Hence by (C3) and the dominated convergence theorem, it follows that 

/ jx(t, x) — t!„(x)j Mo(x)dx —>■ 0 as n —>• oo. 

Jr'^ 

Note that the Theorem of [24] requires the initial probability measure to have finite 
second moments which is trivially met for the auxiliary problems above. 

From all of the estimates above, we can the conclude that 


^cfiYt) 



< \Ecf{Yt)-E^{Xn)\ + 




0 


as n —>■ oo. 


and hence that 

E<j){Yt)= f (f){x)u(t,x)dx for all (^gCj(R‘^). 
JR‘1 
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Now we assume (f> G C'°(R'^) and conclude the proof by a regularization argument. 
Let pe G e > 0, be the standard mollifier and </>£ := pe * ({>■ Then by e.g. 
Section 7.2 in [18], G ||</>e||oo — ll<^llooi 4>e ^ 4> a.e. Observe that 


/R-i 


(j){x)u{t, x)dx — E (^(Yt) 


< 


■>U - / 


+ 


4>eU -E(l)e{Yt) 


+ \EMYt)-Ecl){Yt)\. 


The second term is zero by our result for </> G C'^(R'^), and the two remaining 
terms tend to zero as e —>■ 0+ by the dominated convergence theorem. The proof 
is complete. □ 


4. Proof of Proposition 2.10 

This section is devoted to the proof of Proposition 2.10. The key idea comes 
from the proof of Lemma 16 in [8], and we extend it to the case with unbounded 
coefficients. 


Proof. Since u G there is i? > 0 such that supp u C i?R(0). The PDF 

after one-step Euler-Maruyama scheme is PrU. By Lemma 2.7, 


J- [PtU — u — tAu) 


1 


= / u{x)e^^^ fc(^, x,t) — 1 — T{ib{x)^ — a(a^)0 ) dx 

JRrf V 2 / 

Without loss of generality, we assume d = 1. Note that , we have 

J- [PrU — u — tAu) 


= r^ / (1-0) / e 


10 

-Dl 


/ ( 1 - 0 ) / iPm,uix){i^tO'^i-Stf.'^rdxd9, 

0 „_n 


m,n—0 


where il^m,n ^ C'c(M) and for all x € M, 

4 

\tpm,n{x)\ < Cm.nY l^^'“(^)| ’ 
fc— 0 

and the positive constant Cm,n depends on 6, cr, their derivatives up to fourth order, 
and R. 
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Now we perform inverse Fourier transform and estimate it in L^-norm 


\PtU — U — T^u|| 

= f f / (1 _ 0) j ^ii^gST{ib(x)i-\a{x)f) 
JwL JR Jo Jr 


■ ^ ^Prn,n{x){i6t^)'^{-9t^'^)'^dxd6d^ 

m,n—0 
r*l 


dy 


< r" 


// 




■ tjjm,n{x){iet^)^{-9t^'^)’^dxdC 

m,n—0 
rl 


dyd9 




m,n—0 ^ 

rl 


dydO 


/o JR 
4 


^iT)z ^ib{y+y/9Tz)V^ri—^a{y+'/9Tz)rd 


R JR 


Y i>rn,n{y + V^z){iJ^r])'^{-'rf'Tdzdr] 

m,n=0 


dyd9. 


After the above changes of variables, we continue the estimate 


\PtU — u — tAm 


< P 


til E + 

Jo Jr JRm,n =0 

(y drlj dz 

1 ^ 

I I I Y '^m,n{y + J^z){\f9T)'^ 

Jo Jr Jr _n 


dyd9 


‘ m,n—0 


■ ( £)”^+2" f g^‘(v+b{y+VWTz)\/¥T)r]^-^a{y+\/WTz)ri^ 

\ JR Jv—z 


dyd9 


I I I Y '>Pm,niy + J^z){\/^y 

Jo JR JR _ _n 


‘ m,n—0 


jjm-\-2n 


2'xa{y + JWtz) 


_ (v + b(y + V0Tz)V6T)'^ 

_g 2a(y+\/WTz) 


dz 


dyd9. 


V — Z 
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Now there are no more oscillatory integrals, and therefore we can continue the 
estimate as follows 

\\PtU — u — tAu\\^ 

4 ^1 


< 


, _n 0 »^]R 


m,n—0 ® 


2-Ka{y + VWrz) 


_ {v-\-h{y-\-\/O'rz)\/0T)'^ 
_g 2a(y+\/FTz) 


dydzdO 


= T^ 


< P 



yj2TTa{x) 


(«+b(x)C97)^ 
g 2a{x) 


dxdzdO 


g 2 “(®) dxdzdO 


_n-'O -/r 


m,n=0 


1 


A/27ra(a;) 


1 + 


z + '/Wrb{x) 
yja{x) 


2(m+2n) ^ 


(x + V^fc(x))^ 

e 2 a(x) da;d0. 


Finally we make a further change of variable z = £±2:^^M ^nd obtain 

V a(iic) 

\\PtU — u — Tx 4 u||^ 

^ C'lr^ E /V ( [ ^ + drrdd 

Vjm v27r ^ / / 

4 

< CaT^ E llV'm.nlli 


m,n=0 

4 


<Cr2 E||£’''m||i- 


A;=0 


That is, 


P-rU — U 


— Au 


< C{b, a, u)t. The proof is complete. 


□ 


Remark 4.1. From the above proof we see that the consistency estimate is uniform 
in the elliptic number a > 0 in condition (C2). 


5. Proof of Theorem 3.1 

To prove Theorem 3.1 we need first to show that operator A and its adjoint satisfy 
certain properties, including being dissipative. 
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Proposition 5.1. Assume (Cl). Then the operators A on D{A) := is 

densely defined, dissipative, and closable (possessing a closed extension) in 

Remark 5.2. In view of Theorem 3.5 and Proposition 3.4 we only need A to be 
closable. Alternatively we could have shown directly that A is closed with respect 
to the operator norm ||u||^ := ||m||;^ + as in Section 2 of [15]. 

Proof of Proposition 5.1. The operator A is densely defined since D{A) is dense in 
and it is dissipative by Proposition 3.4. Hence it is closable by Proposition 
II.3.14 in [14]. □ 

On the dual space of the adjoint operator of A with 

domain D{A) = C'“(R'^) is given by the definition below. 

Definition 5.3. The adjoint of A is an operator A* : D{A*) —>■ defined as 

(5.1) f uA* f = j fAu for all / G D{A*) and u G D{A), 

JS.<‘ 

where D{A*) = {/ G L°° :3g G L°° such that / fAu = J gu for all u G D{A)}. 

Remark 5.4. By integration by parts in (2.5), A*f is given by (2.2) for any / G 
Cf (R'^). Moreover / G D{A*) if and only if there exists a. g G such that 

(5.2) A*f = g in 2?' (in distributions). 

Proposition 5.5. Assume (Cl) and (C2). Then the adjoint operator {A*, D{A*)) 
is dissipative on (L°°(]R'^), jj-jj^). 

The proof is quite technical, so we postpone it to the end of this section. Now we 
can employ the following result to prove that A generates a contraction semigroup. 

Proof of Theorem 3.1. Note that A with D{A) = C'))°(R'^) is a bounded densely 
defined operator on L^(IR‘^), and that both A and its adjoint A* are dissipative by 
Propositions 5.1 and 5.5. Hence the closure A generates a contraction semigroup 
on (L^(R.‘^), Ij-Jl]^) by Corollary H.3.17 in [14]. The proof is complete. □ 


In the rest of this section we prove Proposition 5.5. 

Lemma 5.6. Assume (Cl) and (C2), then 

D{A*) C J / G C\R‘^) n L^{R^) :V^f G f] LL(R'^) 

I, l<p<oo 

Remark 5.7. Under different assumptions, similar characterisations of the maximal 
domain can be found in Chapter 3 of [32] and Section 3 of [15]. 


Proof. For any / G D{A*), there exists g G such that (5.2) holds by 

Remark 5.4. By the regularity of the coefficients (see (Cl)), we rewrite (5.1) as 



d ^ d \ d 

f “ 2 ^ \u + fJ2 

i —1 J ^=1 


1 

2 


d \ 

U j 

7 = 1 / 
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On the right hand side, the coefficients of u and diU belong to for all 

1 < p < oo. We can therefore apply Theorem 1.5 in [42] to show that / £ 
for all 1 < p < oo. With this extra regularity and integration by parts, 



On the right hand side, the coefficients of u belong to for all 1 < p < oo 

while the coefficient of diU is zero. By Proposition 1.1 in [42] we then have that 
/ G Wjo’c for all 1 < P < oo. 

Finally, by taking p > d in the Sobolev embedding theorem, see e.g. Theorem 7.10 
in [18], we find that V/ G C'(K‘^). The proof is complete. □ 


The next lemma is crucial. 

Lemma 5.8. Assume f G D{A*) and m := 
sequence {xn} C such that 


sup / < oo. 


Then there exists a 


\x P 

lim/(x„)=m, lim(l + jx„l) lV/(x„)l = 0, and lim ■!——= 0. 

n n n Ti 

Moreover, f(xn) exists and 

V2/(x„) < -Id, 

n 

where Id is the identity matrix on R"^. 


Since the second derivatives are only defined almost everywhere, we need the so- 
called Bony maximum principle [6, 30] to prove the above lemma. We quote the 
version given in Proposition 1.2.12 in [17]. 

Lemma 5.9. Suppose p > d and a function w G achieves a local maxi¬ 

mum at point Xo, then for all ^ G R"^, 

liin essinf V < 0, 

r-s-O \ J 

where Br(xo) denotes the ball of radius r centered at xq. 

The key idea of the following proof now comes from Lemma 2.3 in [23]. 


Proof of Lemma 5.8. For each n G N, we define 

gr^ix) := fix) - 

n 

Hence is bounded from above and tends to — oo when \x\ —>• oo. Therefore there 
exists a;[j for each p„ such that 

gnix'n) = maxp„ =: m„. 

R'i 

Now for arbitrary e > 0, we can find first x^ G R'^ and then TV G N such that 


fixf) > m — e and 


n 


< e for all n > N. 
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Hence for all n > 


m>mn = f{x'^) 



n 


= mBiXgn 


> gn{Xe) = fiXe) - 

> m — 2e. 


I ' 1 ^ 

Therefore m„ —>■ m as n —>■ oo, and hence f{x'^) —>■ m and - > 0 as n ^ oo. 

Next we analyse the derivatives. By definitions of x'^ and gn, 

0 = = V/(4) - ^xl 

Hence, as n —?► oo, 

(1 + KI) |V/K)| = (1 + \x'j)^K\ < i (1 + 3|x;n ^ 0. 


Finally for the second derivatives, 

V^g„(a;) = V'^f{x) - -Id, a.e. 

n 

Since x'„ is maximum point for g„, by Proposition 5.9, 

lim (ess inf V^o„ | < 0, and hence lim (ess inf | < —Id- 

r^0\Br{x'J ^ ) n 

By the definition of the infimum and regularity of /, we can find another sequence 
{xn} such that f{xn) is defined, \xn — x'^\ < 1/n, and 

vV(^«) < -Id- 

n 

Using that / € we can take a further subsequence such that also 

\fiXn) - fix'Jl + |V/(X„) - V/K)| < -. 

n 


Combining all the above estimates, we can conclude that 

lim/(x„) = TO, lim(l + |x„|) |V/(a;„)| = 0, lim-^^^ = 0, 

n n n Tl 

and V^f{xn) < ^Id- The proof is complete. □ 

Proof of Proposition 5.5- Let / G D(A*). We may assume that to := \\f\\^ = 
supRd /. The case to = — inf^d / follows in a similar way and is omitted. 

Note that by (2.2), aij{x) = Yl'k=i 'Xik{x)cFjk and a = (cri, • • • ,On), and then 

d 1 ^ 

A* fix) = bi{x)dj{x) ix)y^fix)akix). 

i=l fe=l 

Let {xn} be the sequence corresponding to / given by Lemma 5.8. From the proof 
it follows that for any fixed A > 0, we can always take 

{xn} C {cc G R'" : \Xfix) - A*fix)\ < \\Xf - A*f\\Y, 
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because the complement of the latter set has zero Lebesgue measure in By 
(Cl) and Lemma 5.8, 

lim/(a;„) = m, \im\b,{xn)dj{xn)\ < limif(l + |x„|) |V/(x„)| = 0, 

n n n 

and 

3 2 3 2 

limsupcTfc (a;„)V^/(x„)(Tfc(x„) < lim - |(Tfe(a;„)| < lim-/C(l + |a;„| ) = 0. 

n n n n n 

Hence, 

limsup(-H*/(x„)) > 0, 

n 

and then 


^ll/lloo =Alim/(x„) < limsup(A/(x„) - A*/(a;„)) < ||(A - A*)/||^ . 

n 

Since / € D{A*) and A > 0 are both arbitrary, A* is a dissipative operator, and 
the proof is complete. □ 


Appendix A. Weak Convergence and Convergence of PDFs 

Now we demonstrate the relationship between L^-convergence of PDFs and weak 
convergence of processes. 

Proposition A.l. Assume and X are random variables with PDFs Un and u 
satisfying 

lim \\un - u|b = 0. 

n—>-oo 

Then A„ converges weakly to X, 

lim |E((.(A„) -E^{X)\ = 0 for all (j) G 


Proof |E(/)(A„) - E^(A)| = l/jj^ <^(a;)(u„ - u){x) dx\ < ||(()||^ ||u„ - uH^. 


□ 


The opposite is not true in general as the following example shows. 

Example A.2. For any n G N, we define 

u„(x) = (1 + sin(2n7rx))x[o,i](x). 

Evidently, > 0 and fg^Vn = 1, and hence is a PDF. 

For any n G N, we take Xn to be a random variable with PDF By the Riemann- 
Lebesgue Lemma, A„ converges weakly to a random variable X with PDF 

= X[o.i](a:^)- 

However, does not converge to v in since 

l|Pra — "clli := / I'Cn — v| = f |sin(2n7rx)| dx = — > 0. 

Jr Jo 
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Appendix B. Proof of Proposition 2.14 


Fix t G [0, T] and define r = t/n for n > \bt\. Without loss of generality, we assume 
that b ^ 0. We start with deriving an expression for the iteration of operators P". 

From the scheme (2.14) and an iteration, we find that 

n —1 

A„ = (l + 6T)"Ao + a^(l + 6T)'^'AP„_fc. 
k=0 

Let 'ipniO s-iid cj}{^) = be the characteristic functions of A„ and an N{0, 1) 

Gaussian variable, respectively. To simplify the notations, let C = C(b, a, T, ug) > 0 
denote various positive constants depending only on b, a, T, and mq, and define 


J)t 


a„:=(l + 6r)- ■.=- l)^ pl := ^ 


a2 2((l + 5r)2"-l) 


2 + 6r 


By the independence of the increments, the definition of AP^, and summation of 
exponents, 

n—1 

'0n(6 = i’oioinO n 4>{0Lk\/T0p) = tpo{otnO(^~ ■ 

k^O 

Note that tpoiO = S'lid ipniO = V^Un{^)- Taking the inverse Fourier 

transform, and after elementary computations using convolution, scaling, and change 
of variables, we find that 


[P^Ug){y)=Un{y)= - * 


1 






{y) 


Ug(x) 


exp 


{y - anxf 


dx. 


Now we estimate the L^-norm of the difference of the exact and the approximate 
densities, \\{Pp — Pt)uo\\]^i. Recall (2.13) and decompose 

P”Mo(y) - Ptuoiy) 


upjx) 

1 1 

/3n ^ 


exp< - 


{y - O^nXf 
r'n 
Ug{x) 


\dx- j 

J Jr 


exp 


Uo{x) 

{y - a„x)^ 


f uojx) / f {y - anxf 

JrJW\ i PI 


Ug{x) 


exp 


exp<; -- 
dx 

{y - anxY 


{y - axf 


dx 


^expj- 


{y - anxY 


Jr x/tt^ 

=: El + p2 + Eg. 

We start with some observations. Let a G 
expansion of ln(l + ax) about x = 0, 

( n / -I . 


- exp < - 


(y - axY 

/32 


dx 

dx 


By I’Hopital’s rule and a Taylor 


1 
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and hence there is a locally bounded function C{a) > 0 such that 



„ / ON" 

~ 

(B.l) 

e“- (1 + -) 

< C(a) — 


V n/ 

n 

Now for El , obviously 

f exol 

^}dx 

yse fdn Since 




( kl 

— 1 < X, 

(B.2) 

|e" - 1| > <^ 1 " 




X < — 1 , 


< ||uo|li- We turn to anal- 


it follows that 


In other words, 




-1 < 2bt, 
2bt < -1. 


(B.3) 


1 ^ IV 
P '/oH 


for all t > 0. 


Moreover 

(B.4) = ^ ^ ^ 

Hence combining (B.4) and (B.l), we get 




P^-Pl 

< 

P^-Pl 

^ C{2bt)AbH'^ Ct 

V/3n p) 


P{P+Pn) 


/32 

~ n 2t n 


Therefore, Hillilli < Ctfn. 

Next, we estimate i? 2 - By a change of variable and the commuting property of 
convolutions, we have 



Oin 

1 




_ e~W 


dx 


— 1 da; 


First we observe with (B.l) and (B.2) that 

q2 -v. 1-1 /I I L^\2n| \ I -I ^2bt\ \ „2bt /-i . t_\2n| 

n>4CT(ivm. 

Hence again from (B.4) and (B.l) we have | ^"^ 2 ^ \ < Therefore 


- 1 < 




- 1 < 


- 1 . 
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So 




-1 


<(1 —e — 1) = e — 1. Nowwe can use Fubini’s 


theorem to conclude that 
1 


mwi < 


Uo 




'-1-0 




-1 


dxdy 


<K||i^ [ -l)d: 

Js. 

= Ikolli ((1 - - 1 


< 


Ct 


for n large enough. 


It remains to analyse E^, and here the computations will differ according to the 
assumptions on uq- We first assume that |a;ito(a^)| dx < oo. Denote 7 „ := ■ 

Combining (B.l) and (B.3), it follows that 

, , ^ ciuye (1V v^) ^ ct3/2 

\ln\ < -^- < 


n y/aH n 

Therefore by a change of variable y y + O'jnX, we have 


l|i?3|li< jj 


Uoix) ( 


V 


exp < - 


(y - anxf 


/32 


- exp < - 


(y - axf 


dxdy 


|uo(a;)| 


^g-(I/+7na:)‘ _ g-y 


dydx 

Jo 


de 


dydx 


^‘Al^ /" \xuD{x)\dx 


< 


Ct3/2 


Izle ^ dz 


for n large enough. 


For the other two cases, we write 


E3= f (—uo(—)--uo(-)) 

JK\an a \aJ) 

1 

ypK^' 


1 1 

Otfi Oi 

1 
a 


J Jr ^OLn' J-K 


1 

e ^ dx 


( ( ^ \ 


Mo (-) 

- Mo ( - ) 

V 7a„/ 

\aJ J 


1 (y-x)^ 

e ^ dx 




1 


=: E'^ + -E'A 

a 


Note that ctn —t a > 0 and hence 0 < a/2 < < 2a for large n. Since lim„ a„ = 

a, and a,an are bounded away from 0. Combining (B.l), it follows that 

^ 1 C{bt)bH^ ^ Cf_ 


1 

1 


e''* - (1 + 6 Vm)” 

a 



OiOCfi 
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Hence ||i? 3 ||i < for large n. 

Regarding to E'^, we assume now \xu'Q(x)\dx < oo. Then similarly 


IS3II1 < 



< 






< C 


1 

Otn 


/ \xur,{x)\dx < -. 

/r n 



dOdx 


Finally we consider the most general case Mq G L^(M) only. We will prove a uniform 
convergence to zero without any rate. For any e > 0, there exist i? > 0 and iVi > 0 
such that 


/ 


l\x\>R 

\aJ 


dx < - and 

O 


'\x\>R 


Uq 


(-) 


dx < - for all n> Ni. 

o 


Next, since translation is continuous in i.e. lim/j_>,o |rc(a;) — w{x + h)\dx = 0 
for any w G and h„{x) := x tends to 0 uniformly for |a:| < R, 

there exists N 2 > 0 such that for n > N 2 , 


[ uo(-)-uo(—)dx=[ uo(-)-uo(- + K{x)) 

J\x\<R ^Oln'' J\x\<R / 

refore w( 

l|i?3 111 < 


dx < 


Therefore we can conclude that E'^ —>■ 0 in L^: 


' |tc|>R 

< e for n > A^i V A^ 2 - 
The proof is complete. 
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/ X \ 
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Wo - 

+ 

Wo - 

o 

3 

1 

1 

o 

3 

\ 

\aJ 


Va„/ 

J J\x\<R 

\a/ 


dx 


References 

1 . L. Angiuli, D. Pallara, and F. Paronetto, Analytic semigroups generated in by second 

order elliptic operators via duality methods, Semigroup Forum 80 (2010), no. 2, 255-271. 

2. D. Applebaum, Levy processes and stochastic calculus, Cambridge Studies in Advanced Math¬ 
ematics, Cambridge University Press, 2004. 

3. V. Bally and L. Caramellino, On the distances between probability density functions, Electron. 
J. Probab. 19 (2014), no. 110, 1-33. 

4. V. Bally and D. Talay, The law of the Euler scheme for stochastic differential equations (II): 
convergence rate of the density, Monte Carlo Methods Appl. 2 (1996), no. 2, 93—128. 

5. H. S. Bhat and R. W. M. A. Madushani, Density tracking by quadrature for stochastic differ¬ 
ential equations, Preprint, https://arxiv.org/abs/1610.09572vl (2016). 

6 . J.-M. Bony, Principe du maximum dans les espaces de Sobolev, Comptes Rendus Acad. Sci. 
Paris, Serie A 265 (1967), 333-336 (French). 

7. N. Bouleau and F. Hirsch, Proprietes d’absolue continuite dans les espaces de dirichlet et 
application aux equations differentielles stochastiques, Seminaire de Probabilites XX 1984/85 
1204 (1986), 131-161 (French). 

8 . Y. A. Butko, M. Grothaus, and O. G. Smolyanov, Feynman formulae and phase space Feyn¬ 
man path integrals for tau-quantization of some Levy-Khintchine type Hamilton functions, J. 
Math. Phys. 57 (2016), no. 2, 023508, 22. 

9. S. Gerrai, Analytic Semigroups and Degenerate Elliptic Operators with Unbounded Coeffi¬ 
cients: A Probabilistic Approach, Journal of Differential Equations 166 (2000), no. 1, 151 - 
174. 






















ON SOLVING SDES BY DENSITY APPROXIMATIONS 


25 


10. W. Chai, A. Naess, and B. J. Leira, Stochastic Dynamic Analysis and Reliability of a Vessel 
Rolling in Random Beam Seas, Journal of Ship Research 59 (2015), no. 2, 113-131. 

11. R. Cont and P. Tankov, Financial modelling with jump processes, Chapman and Hall/CRC 
Financial Mathematics Series, CRC Press, 2003. 

12. M. G. Crandall, H. Ishii, and P.-L. Lions, User’s guide to viscosity solutions of second or¬ 
der partial differential equations, Bull. Amer. Math. Soc. (N.S.) 27 (1992), no. 1, 1-67. 
MR 1118699 

13. G. Da Prato and B. Goldys, Elliptic Operators on Rd with Unbounded Coefficients, Journal 
of Differential Equations 172 (2001), no. 2, 333 - 358. 

14. K. J. Engel, R. Nagel, et ah, One-parameter semigroups for linear evolution equations, Grad¬ 
uate Texts in Mathematics, Springer New York, 2006. 

15. S. Fornaro and L. Lorenzi, Generation results for elliptic operators with unbounded diffusion 
coefficients in Lp- and Cb-spaces, Discrete and Gontinuous Dynamical Systems 18 (2007), 
no. 4, 747-772. 

16. N. Fournier and J. Printems, Absolute continuity for some one-dimensional processes, 
Bernoulli 16 (2010), no. 2, 343-360. 

17. M. G. Garroni and J. L. Menaldi, Second order elliptic integro-differential problems, Chapman 
& Hall/CRC Research Notes in Mathematics Series, CRC Press, 2002. 

18. D. Gilbarg and N. S. Trudinger, Elliptic partial differential equations of second order, Classics 
in Mathematics, Springer Berlin Heidelberg, 2015. 

19. E. Gobet and C. Labart, Sharp estimates for the convergence of the density of the Euler 
scheme in small time, Electronic Communications in Probability 13 (2008), 352—363. 

20. J. Guyon, Euler scheme and tempered distributions, Stoch. Process, their Appl. 116 (2006), 
no. 6, 877-904. 

21. M. Hayashi, A. Kohatsu-Higa, and G. Yuki, Local Holder Continuity Property of the Densities 
of Solutions of SDEs with Singular Coefficients, Journal of Theoretical Probability 26 (2012), 
no. 4, 1117-1134. 

22. M. Hayashi and G. Kohatsu-Higa, A. Yuki, Holder continuity property of the densities of 
SDEs with singular drift coefficients, Electron. J. Probab. 19 (2014), no. 76, 1-22. 

23. E. R. Jakobsen and K. H. Karlsen, Continuous dependence estimates for viscosity solutions 
of fully nonlinear degenerate elliptic equations, Electron. J. Differ. Equations 2002 (2002), 
no. 39, 1-10. 

24. P. E. Kloeden and E. Platen, Numerical Solution of Stochastic Differential Equations, Sto¬ 
chastic Modelling and Applied Probability, Springer Berlin Heidelberg, 1999. 

25. V. Konakov and E. Mammen, Local limit theorems for transition densities of Markov chains 
converging to diffusions, Probability Theory and Related Fields 117 (2000), no. 4, 551-587. 

26. V. Konakov and S. Menozzi, Weak Error for Stable Driven Stochastic Differential Equations: 
Expansion of the Densities, Journal of Theoretical Probability 24 (2010), no. 2, 454—478. 

27. S. Kusuoka and D. Stroock, Applications of the Malliavin Calculus, Part H, J. Fac. Sci. Univ. 
Tokyo, Sect. lA, Math. 32 (1985), 1-76. 

28. A. Lasota and M. G. Mackey, Chaos, fractals, and noise: Stochastic aspects of dynamics. 
Applied Mathematical Sciences, Springer New York, 1994. 

29. V. Linetsky, The Path Integral Approach to Financial Modeling and Options Pricing, Com¬ 
putational Economics 11 (1998), no. 1-2, 129-63. 

30. P.-L. Lions, A remark on Bony maximum principle, Proc. Amer. Math. Soc. 88 (1983), 503- 
508. 

31. C. Litterer and T. Lyons, High order recombination and an application to cubature on Wiener 
space, Ann. Appl. Probab. 22 (2012), no. 4, 1301-1327. MR 2985162 

32. A. Lunardi, Analytic semigroups and optimal regularity in parabolic problems. Modern 
Birkhauser Classics Series, Springer, 2012. 

33. S. Menozzi and V. Lemaire, On Some non Asymptotic Bounds for the Euler Scheme, Electron. 
J. Probab. 15 (2010), no. 53, 1645-1681. 

34. R. Metzler and J. Klafter, The restaurant at the end of the random walk: recent developments 
in the description of anomalous transport by fractional dynamics, J. Phys. A 37 (2004), no. 31, 
R161-R208. 

35. E. Mo, Nonlinear Stochastic Dynamics and Chaos by Numerical Path Integration, Ph.D. 
thesis, Norwegian University of Science and Technology, 2008. 


26 


L. CHEN, E. R. JAKOBSEN, AND A. NAESS 


36. A. Naess and J. M. Johnsen, Response statistics of nonlinear, compliant offshore structures 
by the path integral solution method, Probabilistic Engineering Mechanics 8 (1993), no. 2, 
91-106. 

37. B. 0ksendal, Stochastic differential equations: An introduction with applications, Universi- 
text, Springer Berlin Heidelberg, 2010. 

38. P. Protter, Stochastic integration and differential equations, Stochastic Modelling and Applied 
Probability, Springer Berlin Heidelberg, 2005. 

39. M. Rosa-Clot and S. Taddei, A path integral approach to derivative security pricing, II: 
Numerical methods, Cond-Mat/9901279, 1999. 

40. C. Skaug and A. Naess, Fast and accurate pricing of discretely monitored barrier options by 
numerical path integration, Computational Economics 30 (2007), no. 2, 143-151. 

41. J. S. Yu and Y. K. Lin, Numerical path integration of a non-homogeneous Markov process. 
International Journal of Non-Linear Mechanics 39 (2004), no. 9, 1493-1500, Themes in Non¬ 
linear Stochastic Dynamics. 

42. W. Zhang and J. Bao, Regularity of Very Weak Solutions for Nonhomogeneous Elliptic Equa¬ 
tion, Communications in Contemporary Mathematics 15 (2013), no. 04, 1350012. 

Norwegian University of Science and Technology, NO-7491, Trondheim, Norway 

E-mail address: linghuac@inath.iitnu.no 

Norwegian University of Science and Technology, NO-7491, Trondheim, Norway 

E-mail address: erj@math.ntnu.no 

Norwegian University of Science and Technology, NO-7491, Trondheim, Norway 

E-mail address: arvidn@math.ntnu.no 


